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usually used Jaynes-Cummings Hamiltonian as an initial rotating wave approximation 
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characterized by a single relaxation rate, which is largest at resonance and whose 
expression can be given in closed analytic form. 
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1. Introduction 

In recent years the spin-boson model pQ [3] has experienced a strong revival, as it is 
well suited to describe dissipative and decoherence effects on the dynamics of a two-level 
system (TLS) or qubit coupled to a bath. Crucial for the effects of the environment 
on the dynamics of the TLS is the shape of the spectral density of the harmonic bath. 
It is common to assume an Ohmic spectral density, which is linear in the continuous 
bath modes. In this work we concentrate on a so-called structured bath, for which the 
spectral density is Ohmic at low frequencies but exhibits a Lorenztian-shaped peak at 
a certain frequency Q. It has been shown in [I] that a spin-boson model with such an 
effective spectral density can be exactly mapped on the model of a TLS which is coupled 
to a single harmonic oscillator (HO) of frequency Q, where the latter feels the influence 
of an Ohmic bath. 

Due to its wide applicability the TLS-oscillator system has been object of intense 
research along the years. So it reflects for example the physics of single atoms with 
a large electric dipole moment coupled to the microwave photons of a cavity [5], or 
quantum dots in photonic crystals [6} [7] . More recently the model has received quite 
some attention in the field of quantum computation, where two-level systems are used 
to implement the two logical states of a qubit. We will especially focus on the solid- 
state implementation of such systems. Here, two prominent realizations of a qubit- 
oscillator system are the Cooper-pair box (CPB) [8J - [11] coupled to a transmission 
line resonantor [T2] - [Ej and the Josephson flux qubit [18] read out by a dc-SQUID 
[T9] - [22]. Inspired by experiments with real atoms interacting with a cavity mode, one 
speaks for the CPB case of circuit quantum electrodynamics, as now the CPB plays the 
role of an artificial atom and the waveguide acts as a cavity. From such a setup one 
expects a huge step towards the realization of a quantum computer, as the transmission 
line resonator can be used to couple qubits together [9j [23], store the information of 
qubits or to provide non-demolition read-out schemes [121 H5]- Concerning the flux 
qubit, the read-out usually happens through a damped dc-SQUID, which is inductively 
coupled to the qubit. However, through the SQUID enviromental noise is transferred 
to the qubit leading to decoherence and dissipation within its dynamics. The effect of 
this noise on the qubit depends very much on the strength g of the coupling between 
qubit and SQUID and one faces a conflicting situation. On the one hand one wants a 
strong coupling for a good read-out resolution. On the other hand the coupling should 
be minimized to keep the negative effects of the environment as small as possible. In 
[241 [25] it has been shown that the qubit-SQUID system can be described by a spin- 
boson model with an effective spectral density G e s(u) exhibting a peak at the plasma 
frequency Q of the SQUID. Applying the above mentioned mapping an equivalent point 
of view is to consider the SQUID as an LC-circuit coupled to the Ohmic bath and model 
it as a harmonic oscillator. A detailed description of a nondestructive read-out scheme 
is e.g. given in [26J. 

The spin-boson model can be formally solved using e.g. real-time path integral methods 
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[HE]- However, in order to get closed- form analytical results, approximations must be 
invoked. A quite common one is the so-called weak coupling approximation (WCA), 
which is perturbative in the bath spectral density [2]. However, it has been shown that 
for strong qubit-HO coupling g and for small detuning 5 = Q — A b , where A b is the 
qubit energy splitting, such an approximation breaks down [27J, as coherent exchange 
processes between TLS and oscillator are disregared. For an unbiased qubit the non- 
interacting blip approximation (NIBA) used in [28] -[31] circumvents this problem as 
it is non-perturbative in the coupling g and therefore takes correctly into account the 
influence of the oscillator on the TLS. Moreover, it allows an analytic treatment of 
the dynamics. However, the NIBA is known to break down for a biased qubit at low 
temperatures [U [2j. Another approach, which treats the system non-perturbatively 
in the bath is the flow-equation renormalization method j32j [33], where the spin- 
boson Hamiltonian is diagonalized using infinitesimal unitary transformations. However, 
whithin this approach analytical solutions are difficult to find. Recently a polaron 
transformation was used by Huang et al to obtain analytically the population dynamics 
and confirm the Shiba's relation for an unbiased TLS |34j . 

In the case in which the qubit and the HO are considered as the central quantum system 
being coupled to an Ohmic bath, the numerical, ab-initio quasiadiabatic propagator 
path-integral (QUAPI) method (35J, [36] is a nice tool as it enables to cover both the 
resonant regime, where the oscillator frequency is close to the qubit energy splitting, and 
the dispersive regime with the oscillator being far detuned from the qubit [27J [291 130] - 
Moreover, it can be applied to a biased as well as to an unbiased TLS and therefore be 
used as a testbed for analytical results. For qubits being operated at the degeneracy 
point, which means an unbiased TLS, very often a rotating wave approximation (RWA) 
is applied [12], which is expected to be valid for small detuning and yields as starting 
point the Jaynes-Cummings Hamiltonian [371 [38] . This model was first used to study 
a two-state atom interacting with a single, close to resonance cavity mode of the 
electromagnetic field and predicts e.g. the repeated revival and collapse of Rabi 
oscillations within the atomic excitation probability. By condsidering the TLS-HO 
system in the representation of displaced HO states, Brito et al were able to truncate 
the infinite Hilbert space of this system without loosing the effects of the HO on the 
TLS dynamics [39|. However, so far none of these works could provide an analytical 
expression for the dynamics of the dissipative qubit being valid for zero as well as 
non-zero detuning and for both a biased and unbiased TLS. In this work an analytic 
expression for the dissipative qubit 's dynamics which includes the effects of a finite 
detuning and of a static bias is derived. Specifically, starting from the qubit-HO 
perspective, the eigenvalues and eigenfunctions of the non-dissipative TLS-HO system 
are found approximately using Van-Vleck perturbation theory up to second order in the 
coupling g. Notice that no rotating wave approximation is required. Dissipation effects 
are then evaluated by solving a Born-Markov master equation for the reduced density 
matrix in the system's eigenbasis. 

The structure of the work is as follows. The dissipative TLS-HO Hamiltonian and the 
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main dynamical quantities are introduced in section [21 Inspired by the work of Goorden 
et al [291 EO] , we demonstrate in section [3] how the eigenstates and eigenenergies of the 
non-dissipative Hamiltonian can be found approximately using Van-Vleck perturbation 
theory [1Q1UT]. In this way we can provide an analytical formula for the non-dissipative 
dynamics whichs takes into account the full Hilbert space of the qubit-HO system. After 
that, we show how for low temperatures (k-&T < hfl, HA h ) this infinite Hilbert space 
can be truncated and discuss the relevant contributions of the HO to the dynamics. 
In section H] the influence of the environment is investigated, by looking at solutions of 
the Bloch-Redfield equations. Specifically, analytical expressions fo the TLS dynamics 
are obtained and compared with numerical solutions. The main physical features of 
the coupled TLS-HO system are discussed in section [51 To illustrate the effects of 
counter-rotating terms in the Hamiltonian of the qubit-HO system, which are neglected 
performing a RWA, we compare in section [6] our calculations to results obtained from 
the Jaynes-Cummings model. 

2. The model 

In this section we introduce the Hamiltonian for a qubit coupled through a harmonic 
oscillator to a thermal bath. Further, a formula for the population difference between 
the qubit 's two logical states is derived. 

2.1. The qubit- oscillator-bath system 

To set up the model we consider the Hamiltonian of a qubit-HO system, 7Yqho ; which 
is coupled to an environmental bath, Hb, by the interaction Hamiltonian T~Cob, so that 
the total Hamiltonian becomes 



the Hamiltonian of the TLS/qubit and the harmonic oscillator, and the interaction term 



The Hamiltonian of the TLS is given in the subspace {|L), |R)}, corresponding to a 
clockwise or counterclockwise current in the superconducting loop of a three-junction 
Josephson qubit or more generally to the qubit 's two logical states. In the case of a 
superconducting flux-qubit, the energy bias e can be tuned by an applied external flux, 
<3>ext, and is zero at the so-called degeneracy point. The tunnelling amplitude is described 
by A . For e 3> A the states |L) and |R) are eigenstates of 7Ytls, whereas at the 
degeneracy point those eigenstates are a symmetric and antisymmetric superposition of 
the two logical states. Further, B and are the annihilation and creation operator for 
the HO with frequency fl, and g characterizes the coupling strength. We also introduce 




(1) 



(2) 



'Hint 



hga z {B^ + B). 



(3) 
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the energy splitting hA b = fi\J s 2 + A 2 , between the groundstate |g) and the excited 
state |e) of the TLS. Using the transformation 

m -[-Bin{e/2) cos (6/2) ) (4) 

with tan 6 = —A /e and — | < 6 < |, we obtain the Hamiltonian of the TLS in the 
this basis: Htls — R T (Q)HtlsR(Q) = — ^j^z- The states |R) and |L) become in the 
energy basis 

| J R) = cos(e/2)|g) + sin(e/2)|e), (5) 
\L) =-sin(0/2)|g) + cos(0/2)|e). (6) 

The Hamiltonian 7^ho is diagonal in the eigenbasis {\j)} with j = 0, . . . , oo being 
the occupation number: Hylo = J2j For the eigenbasis of the combined 

Hamiltonian H = 7Ytls + ^ho we write 

{|j>®|g>;|j>®|e» = {|jg>;|je>}. (7) 

Following Caldeira and Leggett [12], we model the environmental influences originating 
from the circuitry surrounding the qubit and the oscillator as a bath of harmonic 
oscillators being coupled bilinearly to the HO. Thus, the environment is described by 
TL-q = Y2k fo^kb^bk and the interaction Hamiltonian is 

Hob = (Bi + B) ]T hv k (b\ + b k ) + (fit + B f ^ h^. (8) 

k k Uk 

The operators b\ and b k are the creation and destruction operator, respectively, for the 
k th bath oscillator, io k is its frequency and v k gives the coupling strength. The whole 
bath can be described by its spectral density, which we consider to be Ohmic: 

k 

In |3] it is shown that the above model is equivalent to that of a TLS being coupled 
directly to a harmonic bath including the single oscillator of frequency Q; i.e., a spin- 
boson model [U [2] with a peaked effective spectral density, 

C 2auQ4 fin) 

eff (ft 2 -u; 2 ) 2 + (2™fi) 2 ' 1 J 

The relation between a and the coupling parameter g between the qubit and the HO is 
g = ft a/ a/ (8k) [231 [25] . This second perspective is suitable for calculating the dynamics 
of the qubit using a path-integral approach, as it was done for example in [31J for the 
case of an unbiased qubit (e = 0). The approach in [HTJ, however, being based on the 
NIBA [2], is not suitable to investigate the low temperature dynamics of a biased TLS. 
Thus, in this manuscript we will consider the TLS and the single oscillator as central 
quantum system and solve the Bloch-Redfield master equations for the density matrix 
of this system, which are valid also for the case of a biased TLS. 
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2.2. The population difference 



The main goal of this work is to determine the dynamics P(t) of the qubit. That means, 
we wish to calculate the population difference 

P(t) = Tr TLS Kp rcd (t)} = (R|p red (t)|R) - (L|p rcd (t)|L) (11) 

between the \R) and \L) states of the qubit. The reduced density matrix of the TLS, 

p rcd (t) = Tr HO {p(t)} = Tr HO Tr B {W{t)} (12) 

is found after tracing out the oscillator and bath degrees of freedom from the total 
density matrix W(t) = e~i nt W(0)eTi' Ht . In turn p(t) = Tr B {W(t)} is the reduced 
density matrix of the qubit-HO system. How to calculate this density matrix will be 



shown later. After some algebra, illustrated in more detail in Appendix A, we arrive at 
an expression for P(t), given in terms of diagonal and off-diagonal elements of p(t) in 
the TLS-HO eigenbasis {|^)}. It reads 



n n,m 
n>m 



where 



Pnn{t) = Y { cos0 (js\n) 2 - (je\n) 2 + 2 sin B (j g\n) (je\n)} p nn (t) 

3 

Pnm{t) = 



(14a) 



cos 6 



(jg\n)(m\jg) - (je\n)(m\je) 



+ sin 9 



(je\n)(m\jg) + (je\m)(n\jg) \^{p nm (t)} 



(146) 



with pnmit) = (n\p(t)\m). How to determine the eigenstates of Hqho is described in 
the next section. 



3. Energy spectrum and dynamics of the non-dissipative TLS-HO system 

In this section we show how to find the eigenvalues of the unperturbed qubit-HO 
Hamiltonian Hqho approximately by using Van-Vleck perturbation theory [lOj HI] ■ The 
idea is to take advantage of the degenerate or doublet structure of the energy spectrum 
of the uncoupled (g = 0) TLS-HO system near resonance, e.g. at A b Q. Then, as 
long as the perturbation is small compared to the energy separation of the different 
doublets, the full Hamiltonian will exhibit a similar spectrum of bundled energy levels. 



3.1. Energy spectrum 

The eigenenergies of the uncoupled TLS-HO system are immediately found by applying 
the Hamiltonian 7io = Htls + T^ho on the eigenstates in (J7J): 

n \jg) = ( + h 3 n) \jg) and H \je)=(^ + h 3 n)\je). (15) 
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The dashed lines in figure [T] show the energy spectrum corresponding to (1151) vs. the 
oscillator frequency fl for the five lowest eigenstates. Except for the groundstate, |0g), 
the states |(j + l)g) and |je) are degenerate in the resonant case (Q = A&). Close to 
resonance the spectrum exhibits a doublet structure. With the coupling being switched 
on, the full Hamiltonian TCquo reads 



n 



QHO 



= R^HquoR = H + H 



hit 



hA b „ 



MIB ] B + hg 



A _ 



[B + B 



(16) 



in the basis {|jg); |je)}. In order to diagonalize the Hamiltonian 7i!qho we consider 7i Int 
as a small perturbation, which is resonable as long as g <C A&, Q. Applying Van-Vleck 
perturbation theory we construct an effective Hamiltonian, 

T^-eS = e l5 7iQHoe~ l5 , (17) 

having the same eigenvalues as 7iIqho but no matrix elements connecting states which 
are far off from degeneracy. Thus, 7i e s will be block-diagonal with all quasi-degenerate 
energy levels being in one common block. As in our case always two states are 
nearly degenerate, each block of 7i e g builds a two-by-two matrix. This can be easily 
diagonalized in order to determine the eigenstates. Following [40, IS] we calculate the 
transformation matrix S up to second order in g. The general formulas for both an 



arbitrary Hamiltonian and 7tIqho are given in Appendix B The only surviving matrix 
elements of the effective Hamiltonian, apart from the ones being of zeroth order in g, 
are 



' je;C?'+l)g 



^eff ) 



hAy/j + 1 with A 



and 



n 



cff 



(2) 
/ je;je 

lleS } 

' ig;ig 



fe 2 



A 2 b n 



9 2 +J- 



hAl 



A 2 b (A b + n) 

hiWx + U + iWo}. 



g 2 = K( Wl - jWo 



(19) 
(20) 



Thus, H c g = 7i5l + H^fJ + HplJ has the matrix structure 



(2) 





( •.. 




\ 






f +jtl + Wt-jW 
y/j + 1A 


VT- 

-f + (j + i)n -\ 


-1A 

- W 1 + (j + 2)W 






V 







(21) 

where the section shown corresponds to the basis states \je) and \(J + l)g). From 
this form it is easy to calculate the eigenstates and eigenenergies. The groundstate 
|0)cfr = |0g), which is an eigenstate of 7i e s, has the eigenenergy 



Ao h [ + W + W 1 



(22) 
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The other eigenstates of 7Y e fr are, j > 0, 



|2j + l) eff = cos ( - J 



a 



L ) I(J + 1)S> 



+ sin 



a, 



\2j + 2) eff = - sin (^) |(j + l)g) + cos 



|je), 
b'e), 



corresponding to the eigenenergies 

E2j+l/2j+2 = 



(j + + w 1 + w q t -t^- 

2 2 cos aj j 



h 



(j + ~)n + W 1 + W t-J5* + 4(j + 1)| A 



with 5j = A b - Q - 2(j + 1)W , tana,,- 



2v7TT|A| 



and 



(23a) 
(235) 



(24) 



< atj < 7T. (25) 



By construction these are also eigenenergies of Hqho Using the transformation ([T] 
we get the eigenvectors of 7iIqho as 

|0) = e- i5 |0) eff , |2j + 1) = e~ i5 |2j + l) cfr and |2j + 2) = e - i5 |2j + 2) eff . (26) 

The energy spectrum of 7Yqho is shown in figure [1] for the case of an unbiased TLS 
(e = 0). We want to emphasize that our findings are also valid for the more general case 
e 7^ 0. At resonance, where the spectrum for the uncoupled case is degenerated, avoided 
crossings can be seen. The gap between two formerly degenerated levels for Q = A& is 

(27) 



E, 



'*j+2 ~ E 2j+ i = 2Vi + 10 + 0(g" 

which is as predicted by the Jaynes-Cumming model [371 EH]- As we will show in 
section [6], the second order correction Wq in (12II) . whichleads to a shift in the resonance 
frequency, is a result of the counter- rotating terms in Hquo- As such it can be 
interpreted as a Bloch-Siegert shift 



3.2. Dynamics of the qubit for the non- dissipative case 

With the coupling to the bath being turned off, the time evolution of the density matrix 
of the qubit-HO system is given by p(t) = e~^ HQHOt p(0)e^ WQHO * and consequently 

Pnm{t) = (n\p(t)\m) = e- iaW pnm(0) 

with uj nm = \{E n — E m ). With that (}T3"]) becomes 



Prim 

(0) cos^ nm t, 



(28) 
(29) 



n>m 



where we defined po = J2 n P nn (®)- F rom fl29l we notice that the dynamics of the qubit 
is characterized by an infinite number of oscillation frequencies rather than showing 
Rabi oscillations with a single distinct frequency. This is clearly a consequence of the 
coupling of the HO to the TLS. Further we assume that at t = the qubit starts in 




Figure 1. Energy spectrum of the coupled TLS-HO system vs. the oscillator frequency 
f2. Solid lines show the energy levels for the five lowest energy states with the coupling 
being switched on (g — 0.18) and the TLS being unbiased (e = 0). Frequencies and 
energies are given in units of Ao and hAo, respectively. For comparison the energy 
levels for the uncoupled case are also given (dashed lines). At resonance (Q = Ab) 
the spectrum exhibits avoided crossings, whereas it approaches the uncoupled case 
away from resonance. The vertical dashed lines visualize three different situations: 
the negatively detuned regime (line a), the resonant case (line b) and the positively 
detuned regime (line c). 



the state \R) and that the occupation numbers of the HO are Boltzmann distributed, 
so that 



1 



p(0) = iRXRPe"^ 



(30) 



where Z = e h P n / 2 / (1 — e /3nn ) is the partition function of the oscillator and (3 = (ksT) 1 
denotes the inverse temperature of the system. In the TLS-HO eigenbasis this becomes 



Pnm(0) = (n\p(0)\m) 



j=0 



< '' (,s I -r I (n\jg) + sin (-- ) (n\je) 



I — I (j'sl m ) + s i n f ~W ) e l m ) 



(31) 



3.3. Low temperature approximation 



With ( 1291) we found a formula which describes using the approximate eigenenergies 
and eigenstates in (T24|) and (1261) the non-dissipative dynamics up to second order in g, 
thereby taking into account all oscillator levels. Thus, we still have to deal with an 
infinite Hilbert space. Typically experiments, see e.g. in [T3l [21] . run in a temperature 
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regime for which < Ml, hA^,. Considering the exponential function in (1311) we 
assume the higher oscillator levels to be only sparsely populated and the maximum 
value of the sum in ( 13TI) is truncated to j = 1. Nevertheless, states | jg/e) with j > 1 
still play a role in the dynamics. In fact, due to the Van-Vleck transformation exp(— iS), 
for example the state 

|8> = e- i5 |8) eff = e~ iS [- sin (^) |4g) + cos (^) |3e>] (32) 

yields nonvanishing contributions to the matrix elements (n\lg) and (n|le) occurring in 
fIBTj) due to the fact that the energy eigenstates (|2"6T) of the coupled TLS-HO system are 
made of linear combinations which involve also these states. 

Using (1146j) together with (|28|) one finds that coefficients p nm (0) with n > 7 are of higher 
than second order in g. The same is valid for p 50 , pm, P55 and pqq. Thus, those terms 
play no role in our calculation of P(t). Furthermore, e~ 1 2 /3n (g/At,Q) 2 <C 1. Neglecting 
also these contributions we find that p n , m "C 1 for n > 5. In the end it will be sufficient 
to concentrate on eigenstates of Hquo up to |4). This trunctation leaves us with ten 
possible oscillation frequencies u) nm , where n, m = 0, 1, . . . , 4 and n> m. 

As an example we calculate the dynamics of an unbiased TLS (e = 0). Here the 
coefficients p , p 3 o(0), P4o(0), £>2i(0) and ^43(0) vanish due to symmetry, so that 

P(t) = P10 cos (u w t) + p 20 cos (u 20 t) + psx cos (uj 31 t) + p 41 COS (Uftt) ^ 

+ p 32 COS (U 32 t) + p 4 2 COS (i04 2 t) . 

Additionally as a benchmark we consider the mostly studied resonant case, where 
Q = Aft = A . In this case we find with ff24|) the transitions frequencies 

u w = A - g, uj 20 = A + g, (34a) 
u 31 = A + (1 - V2)g, u A1 = A + (1 + y/2)g, (346) 
u 32 = A - (1 + y/2)g, u 42 = A - (1 - V2)g. (34c) 

The dynamical quantity P(t) and its Fourier transform are shown in figure El One 
clearly sees the influence of the coupling to the HO on the dynamics of the TLS. 
Instead of Rabi oscillations with a single frequency, P(t) oscillates with six different 
frequencies, which are in the Fourier spectrum symmetrically located around the point 
ijj = A . Among those frequencies c^io and u 20 are dominating. They correspond 
to transitions between the first or second energy level of the qubit-HO system and its 
groundstate and their weight is almost equal. To summarize, one notices that due to the 
coupling with the oscillator additional frequencies are induced into the qubit dynamics. 
Theoretically, the number of those frequencies is infinite. At low temperatures, however, 
transitions between the lower energy levels of the system are clearly dominating. Again, 
for simplicity we have shown here the case of an unbiased TLS. For e 7^ the behaviour 
is similar only that in the Fourier spectrum the weight difference of the two dominating 
peaks will be more pronounced. 



Dissipative dynamics of a biased qubit coupled to a harmonic oscillator 



11 




Figure 2. Left-hand graph: Dynamics of the population difference P(t) for the 
unbiased TLS-HO system at resonance (SI = Ao) with g = 0.18 and j3 — 10. Right- 
hand graph: Fourier transform F(u) of P(t) for the unbiased system. The peaks are 
situated around u> = Ao according to (|34a[) - (134 c[) . Clearly, u>io and U20 are the 
dominating frequencies. In order to visualize the delta- functions, finite widths have 
been artificially introduced. 



4. The influence of the environment 



In the preceding section we neglected the influence of the bath on the qubit-HO system. 
Yet, in order to model a realistic situation, we have to pay attention to environmental 
influences, as they lead to decoherence and dissipation in the dynamics of the qubit, 
which is harmful for quantum computing application. Thus, we will now consider the 
full Hamiltonian 7i. 



4-1. Master equation for the qubit-HO system 



As shown in section 12.21 we need for the calculation of the population difference P{t) 
the density matrix p(t) of the qubit-HO system. Starting from the Liouville equation of 
motion for the full density matrix W(t) of 7i, 

ih^^ = [H B(t) h W(t) l ], (35) 

where the index stands for the interaction picture and following [121111], we can provide 
a Born-Markov master equation for p{t) being in the Schrodinger picture and expressed 
in the basis of the eigenstates of 7Yqho : 

win. klPk,l (t). (36) 

kl 

The free dynamics of the system is given by the first term of the right-hand side in the 
above equation. The rate coefficients are defined as 

£-nm,kl = [G(tO n k)N n k — G{UJl m )N m l] X n kXi m 

—8mi G(uJi'k)NirkX n i>Xi>k + 5 n k G(iOik')Nk>iXik'Xk'm ^ 
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with N nm = |[coth(^,/3o; nm /2) — 1] and X nm = (n\ (B^ + B) \m). For the derivation 
of the master equation, besides the Born and Markov approximations, some more 
assumptions have been made, which we briefly mention. First, we consider our system 
and the bath to be initially (at t = 0) uncorrelated; i.e., Vl^(0)i = p(0)pb(0) with 
Pb(0) = Zg 1 e _/3WB and Z-q the partition function of the bath. Further, with the bath 
consisting of infinite degrees of freedom, we assume the effects of the interaction with 
the qubit-HO system to dissipate away quickly, such that the bath remains in thermal 
equilibrium for all times t: W(t)\ = p(t)ipB(0). Besides, an initial slip term which 
occurs due to the sudden coupling of the system to the bath is neglected [2]. And as 
last approximation the Lamb-shift of the oscillation frequencies u nm was not taken into 
account [H]. 



4-2. Matrix elements 

In ( 1371) X nm describes matrix elements of the operator X = (B + B^) in the qubit- 
HO eigenbasis. By use of (123 ap . (12351) and (I26I) those states were expressed in the 
basis {|jg); |je)}, and we will also calculate the oscillator matrix elements in this 
basis. For that purpose the operator X = e lS (i?^ + B) e~ lS is defined. Four different 
situations can be distinguished. There are matrix elements were neither the oscillator 
nor the qubit state are changed, namely (jg\X\jg) = — 2L and (je\X\je) = 2L Q 
with L = eg/A h Q. We see that those elements are independent of j, the occupation 
number of the oscillator. Next, we look at the case where a single quantum is emitted 
or absorbed from the oscillator and get (jg\X\(j + l)g) = ^/j + 1(1 + L osc ) and 
(je\X\{j + l)e) = VJTT(1 - Lose) with 

(2A 6 + 3Q)Ag 

Losc ~ AiQ(A b + ny 9 ■ {6 * } 

For a transition within the qubit we have (jg|X|je) = A g / A^Ab + Q) = L q . And 
finally, if the qubit and the oscillator state are changed simultaneously, one obtains 
(jg\X\(j + l)e) = v/jTTL+ osc and (je\X\(j + l)g) = v^TTL q>osc , 

where L —= At(A b + m^ + m 9 and L — = Atn(A b -2n) 9 - (39) 

Comparing the magnitude of the transition terms, we notice that those consisting in 
changes of the oscillator occupation only are the dominant ones, as they have a part 
which is of zeroth order in g. Further, for the case in which the qubit is operated at the 
degeneracy point L and Lctosc vanish. With those results we can calculate the matrix 



elements X nm . They are given in Appendix C 



4-3. Dynamics in the dissipative case 

Like in section 13.31 we assume the system to be operated at low temperatures and thus 
take as highest qubit-HO state the eigenstate |4). For determination of P(t) the formulas 
of section 12.21 can be used. Unlike in the non-dissipative case pit) is not given anymore 
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by the simple expression (1281) . Rather we have to solve a system of coupled differential 
equations, namely fl36l) . To do this analytically we will follow three different approaches 
and compare them finally to the numerical solution of ( 1361) . We start by introducing 

Pnmit) = e- iW -V nm (t), (40) 

which yields the set of differential equations for a nm : 

0-nm{t) = *J2 ZnmM^^^Mt)- (41) 
kl 

4-3.1. Full secular approximation (FSA): As a first approach we make the full secular 
approximation; i.e., we neglect fast rotating terms in (T4T]) and keep only contributions 
where cu nm — u> k i vanishes. In this way the off-diagonal elements of a nm are decoupled 
from the diagonal ones so that 

0~nn(t) =7l ^2 £nn,kk a kk(t) , (42) 
k 

<7nm{t) = KC nm>nm a nm {t) for n^m. (43) 
The equation for the off-diagonal elements is then 

° nm {t) = a° nm e^—\ (44) 
which becomes with ([40]) 

n 1 4-\ n „irCnm,nmt e — lU nrn t { A^} 

Fnm\ L ) — Hnm c e ■ \^°) 

As through the FSA the oscillatory motion of the dynamics is separated from the 
relaxation one we can divide (|T3|) into two parts, 

P(t) =Prelax.(t)+ Cephas. (*), (46) 

where P rc iax.(^) = J2 n Pnn{t) describes the relaxation and P d e P ha S .(t) = E„> m Pnm(*) tne 
dephasing parts of the dynamics. With fj45|) the latter takes the form 

Pdcpha S .(t) = Yl Pnm(Q) e ~ rnmt COs(u nm t) (47) 

n>m 

with the dephasing rates T nm = — 7r£„ m nm . Expressions for the dephasing coefficients 
£>nm,nm can be found in Appendix D and the initial conditions p° m = cr° m = p nm (0) 
are given by fl3TT) . The diagonal elements are more difficult to obtain, as one has to 
solve a system of coupled differential equations, ( 1421 . Calculating the corresponding 
rate coefficients of this system for the five lowest eigenstates, we find that there are 
only eight independent ones, namely £oo,u> £00,22, £11,22, £11,33, £u,44, £22,33, £22,44 
and £33,44. They are given by 

C jjM = 2G(u jk )N jk Xj k with j < k, (48) 

where j and k adopt the above values. Furthermore, £00,33, £00,44, £33,00 and £44,00 
vanish. The remaining rate coefficients are combinations of the above. We find that 

£fcfc,jj = £/j,fcfc + ^(uj k )Xj k = (Njk + l)2G(u}j k )Xj k (49) 
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and 



(50a) 
(506) 
(50 c) 
(50 d) 
(50e) 

However, the system (|42p is still too complicated to be solved analytically. Thus, we 
invoke a further approximation: we consider the factor N nm + 1 = ^[coth.(h/3u nm /2) + 1] 



£00,00 — 


— £ll,00 


— £22,00 






£ll,ll = 


— £oO,ll 


— £22,11 


— £33,11 


— £44,11 


£22,22 — 


— £oO,22 


— £ll,22 


~~ £33,22 


— £44,22 


£33,33 = 


— £ll, 33 


— £22,33 


— £44,33- 




£44,44 = 


— ^11,44 


— £22,44 


— £33,44- 





with n < m in (09]) and use that lim^^oo coth(ft/3u; nm /2) = — 1. It depends strongly on 
the temperature j3 for which value of u nm this limit is reached approximately. For the 
parameters we are working with one usually is in the region where (N nm + 1) <C 1. Thus, 
we will neglect in the following terms containing the factor (N nm + 1). Furthermore, one 
sees from ( 1271) that uj\2 ^ 9 and CJ34 ^ g. With that £11,22 = 0(g 3 ) and £33,44 = 0(g 3 ) 
can be neglected. Using (150 a\i - (150 eft the matrix of the system (H2|) becomes 



£ 



relax. 



/o 


£00,11 


£00,22 








\ 





— £00,11 





£11,33 


£11,44 










— £00,22 


£22,33 


£22,44 













— £11,33 — £22,33 







\o 











—£11,44 — £22,44 


/ 



(51) 



The eigenvalues and eigenvectors of this matrix and the associated time evolution of 



the elements a nn (t) are given in ( 15. ip - (15. 5p of Appendix E Unlike for the dephasing 
part (|4Tp . we cannot extract a simple analytical expression for the relaxation rate as 
Preiax.it) = J2 n P"ri(t) now consists of a sum of several exponential functions, cf. ( I14ap 
together with ( 15. ip - (15.51) . But still we are able to provide an analytical formula for 
P(t) using (US]). 



4-3.2. An ansatz for the long-time dynamics: In order to obtain a simple expression for 
the relaxation part, we consider the long-time dynamics of the system. In other words, 
rather than looking at the many relaxation contributions to the populations o~ nn (t), we 
focus on the smallest eigenvalue of the relaxation coefficients, as it will dominate at 
long times. Further, we consider only the rate matrix associated to the three lowest 
qubit-HO eigenstates, |0), |1) and |2) in ( j42l and obtain with ( 150aP - ( 150 d that 

^ — £11,00 — £22,00 £00,11 £00,22 » 

£rclax. 



V 



"£11,00 ~ £22,00 
£11,00 
£22,00 



£00,11 
—£00,11 — £22,11 
£22,11 



£00,22 
£11,22 
— £00,22 — £11,22 



(52) 



Here, we have not neglected the rate coefficients containing the term (N nm + 1) and 
further took £11,22 into account despite of being of third order in g as such contribution 
removes the degeneracy between the two lowest eigenvalues at resonance, cf. inset in 
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nn,mm 



15 



^ 4(A)0, ll^OO, 22 + ^11,00^00,22 

+^00,11^11,22 + ^11,00^11,22 + ^00,11^22,00 + ^11,22-^22,00 

l 1/2 \ 

+^22,11^00,22 + ^11,00^22,11 + ^22,00^22,ll) f • 

With the system being detuned this expression can be simplified further, namely 



r ?« 

J. r ~ 



7lC 



00,22 



for tt < A b ; T r « tt£ 



00,11 



for Q > A b . 



(53) 
(54) 



0,15 



3" 




0,05 



Figure 3. The relaxtion rate I\ as it is given in (|53|) drawn against the oscillator 
frequency f2 (solid line). Used values are e = 0.5Ao, corresponding to a frequency 
splitting Ab = 1.118A , and coupling g = 0.18A . Moreover, the damping constant 
is k — 0.0154 and (3 — 10(fi.Ao) _1 . At resonance (ft = Ab) T r is maximal. For a 
comparison also the second smallest eigenvalue is plotted (dashed line). The inset 
shows the two eigenvalues close to resonance. 

In figure [3] the relaxation rate r r as obtained from ( J53l) is shown as a function of the 
oscillator frequency Q. Clearly, it is maximal at resonance (Q = A b ), whereas it decays 
with Q being detuned from the resonance. This effect has already been predicted by 
Blais et al [12J. As the qubit is not directly coupled to the bath but rather through 
the oscillator, the latter being detuned filters out the environmental noise at the qubit 
transition frequency. Additionally, we show the second smallest eigenvalue of (1521) . We 
notice that close to the resonant point (Q = A b ) there is an avoided crossing. Finally, 
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P 



we find that 

-Prelax.(^) = (PO ~ Poo)e~ Vrt + Poo, (55) 

where like in section l3T2"l p n = ^ n Pnn(0). For getting p^ we have in principle to find the 
steady-state solution of (1121) . Here, we just assume for t —>■ oo a Boltzmann distribution 
for the qubit-HO system, so that p nn (oo) = Z^^e -13 ^ with -Zqho = Sn e_/3£; ™- Thus, 

^2^2 {cos B (jg\n) 2 - (je\n) 2 + 2sin9(jg|ri)(je|n)} p„„(oo). (56) 

n i 

The formula for the long-time dynamics is obtained as 

P(t) = (p - Poo)e~ rr * + Poo + Yl Pnm(0)e' Vnmt cos(uj nm t). (57) 

n>m 

To get further insight on the dominant frequencies we evaluate the Fourier transform of 
( [571) according to 

POO 

F(u) = 2 dt cos utP(t), (58) 
Jo 

yielding 

r r 

F(u) = 2(p - Poo) 2 ; ^ 2 + 27rp 00 <S(u;) 

CJ + 1 r 



1 

+ 



L^L + {Umn + U) 2 + [UJ mn — W) 2 



(59) 



4-3.3. Partial secular approximation (PSA): An improvement to the FSA is to take 
into account certain non- vanishing contributions of to n m~ ljJ u- We have to keep in mind, 
that there are quasi-degenerate levels close to resonance. In our case the first with 
second energy level and the third with fourth one build a doublet, meaning that they 
are close together in energy space. The level spacing is approximately proportional to 
g for the former and \/2g for the latter. Because of that and as the transitions from 
level three and four are less probable, we will in the following only consider the first and 
second level as being almost degenerate. Taking this into account in (j4ip we arrive for 
the diagonal elements at 

0-nn(t) =^Y1 C nn,kkO- k k{t) + vCnn^U (t)e"^ 12 ' + 7rC nn>2 l(T21 (t)e -iW!llt . (60) 
k 

A numerical analysis shows that the effect of the last two terms on the right-hand side 
of the above equation will in the worst case lead to very small wiggles in a nn (t) and play 
no role in P(t). Thus, we finally write 

kkO~kk{t), (61) 

k 

which is the same equation as we got in the FSA approach. However, the off-diagonal 
contributions croi, 002, cr 13 , <7 2 3, au and 024 have to be examined more carefully. From 
(j4Tj) we find that one has to solve the equations 

Pnrait) ( i^Vim ~\~ Cnm,nm) Pnm(t) ~\~ 7T C-nm,jkPjk{t) , (62) 
Pjk{t) = n£jk,nmPnm{t) + (lUJjk + 7T £jk,jk) Pjk{t) (63) 
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with {(nm), (Jk)} = {(01); (02)}, {(13); (23)} or {(14); (24)}. As solution one gets 



(+) „.(+) „A<t> * , „H ..(-) „A 



(-) 



r\ = ' ny- 1 p nm,jk _J_ ' <ii* ' p nm,jk 

Hum ^nm,jk nmjk 1 nmjk nmjk i 



,(+) 



(+) 



(-) 



(64) 

Here, the oscillation frequencies and the decay of the off-diagonal elements are given by 



with 



Rnmjk Y ["'(Aim,ljm £jk,jk) '(^nm ^jfc)] "I - 47T £>nm,jk£jk,nm' 

The amplitudes of the oscillations are given through the coefficients 

27T £jk,nmPriTn Pjk ["'(^nm 1 Tim £jk,jk) '(^nm ^jk) ~F ^nrajfc] 

2-Rnm j'fc 



J+/-) 
nm,jk 

and 

(+) 



2tt 



nm,jk £ 



■jk,nm 



ram ^jfc 



nm.jkj 



(66) 
(67) 

(68) 
(69) 




0,6 0. 



0,075 



0,05 



0,025 



Figure 4. Comparison between the dephasing rates of the two dominant frequencies as 
they are obtained using the FSA or the PSA, respectively. On the left: the FSA rates 
Toi = — 7r£oi,oi (dashed red line) and To2 = — 7r£o2,02 (solid black line). On the right: 
the real part of Aq7 02 ( re ^ dashed line) and Aqi 02 (black solid line) as given by ([55)1 is 



shown. The rate dominating the dephasing behaviour is defined as = ^t{Xoi' 02 }- 
For fl < Ao we see that i s approximated by the FSA rate T02, while for fl > Aq 
by r i. Used values are e = 0, g = 0.18, K = 0.0154 and (5 = lO(ftAo)- 1 . 



Thus, we have again all ingredients to calculate analytically the relaxation and dephasing 
part of f)46p . For the PSA we cannot provide a simple expression for the dephasing 
rates as in the FSA, where we had Y nm = — nC nm ^ nm . As one can see from (1641) and 
f l65|) . poi and p 2 are a mixture of contributions decaying with 3?{Aq^o 2 } and 3?{A 7o 2 }- 
Similar to our findings for the relaxation rate, also here the smallest eigenvalue will 
dominate the dephasing behaviour. From the right graph in figure H] we find that 



this is 3?{Aqi 2} = r^2 j . Comparing it with the dephasing rates we got using the 



,(+) 
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FSA, left graph in figure IU we see that for negative detuning (fl < A b ) the rate 
To2 = _ 7r ^02,02 approximates T\t,\ whereas for positive detuning (Q > A b ) this is 
done by r i = — 7r£ i,oi- In the FSA T 2 and r i correspond to the frequencies u w and 
co>20) respectively. In the PSA the frequency id^ = S{Aq^q 2 } is given by oo[p = oj 2 o 
for < A b and u^' = co> 10 for > A b . Hence, for negative detuning oscillations with 
frequency u 2 q will dominate the dynamics, while those with u>io will almost vanish. For 
positive detuning it is the other way round. In (1M1) and (1651) this behaviour is reflected 
by the coefficients c^mjk anc ^ v nmjk- Around resonance (Q ~ A b ) the PSA tells us by 
( 164|) and (1651) that the dephasing rates and frequencies are a mixture of r i and r 2 or 
uio and u 2 o, respectively. From the left graph in figure H] one notices further that the 
FSA rate To2 grows linearly with Q for positive detuning. However, as the weight of the 
corresponding frequency u 20 will be almost zero, r 02 will give no relevant contribution 
to Pdephas.it) in this regime but the dephasing will rather be associated to the FSA rate 
r i- Hence, out of resonance the FSA will still fairly well describe the dynamics of 
P{t). Comparing the expressions for £01,01 and £02,02 given in 



Appendix D by ( 14.ip and 



( 14. 2p with the approximative expressions for the relaxation rate at positive and negative 
detuning (154]) . we see that for zero bias (e = 0) the PSA dephasing rate is equal to T T /2. 
For a biased system an additional term is added depending on the spectral density of 
the bath at uj = 0. In figure [5] we compare the three analytical solutions described 




P(t) numerical 
P(t) FSA 
P(t) PSA 
P(t) Ansatz 



100 
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F(co) Ansatz 
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Figure 5. Comparison between the behaviour of P(t) and its Fourier transform F(w) 
as obtained from the numerically exact solution (black solid curve) of the equation (l36|) 
and the three analytical approximations discussed in the text. The red dashed curve 
is the full secular approximation (FSA) solution, the green dotted-dashed curve the 
partial secular approximation (PSA) solution and the blue double-dotted-dashed curve 
the analytical formulas (|57|) and (|59| . The parameters are e = 0, Q = Ao, g = O.I8A0, 
K = 0.0154 and j3 = lO(fiAo)" 1 . For the chooscn regime of parameters differences 
between numerical and analytical results are barely visible. 

above to the numerical solution of the master equation for the case of an unbiased TLS 
being at resonance with the oscillator. Concerning both the dynamics of P(t) and its 
Fourier spectrum we see a good agreement between the different solutions. The one 
being closest to the numerical solution is the PSA solution. We also want to mention 
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that going to stronger damping k, the FSA results start to show deviations from the 
numerical solution. Here, one should use the PSA only. However, for the parameter 
regime used in the following, we will mainly apply ( 1571) due to its simple, analytical 
form. 



5. Discussion of the results 



Having solved the master equation fl36l) analytically and numerically we can examine 
the dynamics of the system and its Fourier transform for different situations. First, 
we will look at a qubit operated at the degeneracy point (e = 0) being in and out of 
resonance with the oscillator. Then, we will concentrate on the biased qubit in the same 
regime of parameters. 



5.1. The unbiased qubit 

For unbiased qubits we can compare our predictions with the analytical results obtained 
in [31] by starting from a spin-boson model with the effective spectral density (|T0|) . In 
[3T] a so-called weak damping approximation (WDA) based on the non-interacting blip 
approximation (NIBA) is applied. The WDA allows a non-perturbative treatment of 
the coupling between the TLS and HO and hence can reproduce the occurence of two 
dominating frequencies as expected e.g. from exact QUAPI calculations [27]. The 
NIBA, and hence the WDA, however, become not reliable for a biased TLS. We find 
that the overall agreement between our approach and the WDA is very good. However, 
in the WDA solution the frequencies are slightly shifted compared to the ones obtained 
from our master equation. This may result from the perturbative expansion we have 
performed with respect to g by applying the Van-Vleck perturbation theory. 
First, we look at the resonant case shown in figure El In agreement with previous works 
[271 13T] . we find that the dynamics is dominated by two frequencies corresponding to 
uj w and CJ20 with separation being approximately 2g. The weight of the latter is a bit 
larger. The reason for the bigger weight is that at resonance (Q = Ab) the qubit- 
HO eigenstate \j) is not a symmetric superposition of the states \j, e) and \j + l,g) 
unlike it is predicted by the Janyes-Cummings model (cf e.g. [12]). We notice that 
the two unequal peaks have indeed been experimentally observed in [13] (see Fig. 4b 
therein). Considering the states |l) e fi and |2) efr in (123 qj) and (123 6j) . one already sees 



that for a symmetric superposition of these states we need that 5q vanishes or that 
Q = [{Al + 2g 2 Al)/Al}-^ (cf . (1751) ). Besides, in order to get the qubit-HO eigenstates 
one still has to perform the Van-Vleck transformation, which adds contributions to 1 1) 
and 1 2) from states corresponding to oscillator levels higher than j = 1. Thus, our 
system behaves for Q = A as being negatively detuned, which means that the peak 
belonging to the higher frequency dominates, as we will show below. Slightly increasing 
Q will give a stronger weight to the peak at uj 10 . This effect is not very pronounced 
for the non-dissipative dynamics of the unbiased qubit (figure [2]), as there the two 
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Figure 6. Dynamics of the population difference P(t) (left-hand side) and its Fourier 
transform F(cu) (right-hand side) for an unbiased TLS being in resonance with the 
oscillator (fi = Ao). The numerical solution of the master equation (black solid line) is 
compared with the analytical expressions ([57]) and (|59[) (red dashed line) and the weak 
damping approximation (WDA) from [31j (green dotted-dashed line) . The parameters 
are e = 0, fl = A , g = 0.18A , k = 0.0154 and (3 = 10(7iAo) -1 . From the Fourier 
transform one sees that two frequencies are dominating the dynamics. The separation 
of those two peaks is approximately 2g. The non-dissipative dynamics is shown for 
comparison in figure [H 



frequencies are still almost equally weighted. Looking however at the Fourier transform 
of the dissipative dynamics (|59|) . one notices that the relaxation rate also contributes 
to the weight of the peaks with a prefactor T~^. As for a negative detuned system r i 
is slightly bigger than r 2, the difference between the two peaks becomes more clear in 
the dissipative case. For e ^ the effect can already be noticed in the non-dissipative 
case. 

Next, we consider in figure [7] the case of negative detuning, where Q < Ao- No matter 
which approach one is looking at, clearly the frequency U20 is dominating. Furthermore, 
paying attention to the timescale of the dynamics, one notices that the relaxation time 
is enhanced compared to the one we found for the resonant system. This behaviour was 
already explained by the formula (153|) for the relaxation rate. Again, the numerical and 
the solution obtained by using the long-time ansatz in section 14.3.21 agree quite well with 
each other, whereas the amplitude of the oscillation with frequency u>2o is stronger in 
the WDA approach. Also remarkable is the fact that looking at the Fourier transform 
in figure [7] one sees in the inset already small contributions of the higher oscillator levels. 
The transitions corresponding to co> 2 4 and co> 2 3 give raise to small additional peaks, while 
the contributions of U13 and Uu are negatively weighted and cause dips. The WDA 
approach does not show this additional contributions. They are, however, confirmed by 
the numerical QUAPI calculations in |27j (see figure 2 therein). In the case of positive 
detuning (fl > A ) shown in figure M we find a quite good agreement between all three 
approaches. Also for postive detuning the relaxation time is enhanced compared to the 
resonant case. Contrary to the negatively detuned situation the additional peaks have 
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Figure 7. Dynamics of P(t) and its Fourier transform F(u>) for negative detuning 
(0 < Ao) and for e = 0. Same parameters as in figure [5] are used except that now 
il = 0.75Ao. The frequency W20 dominates the dynamics. The inset on the right graph 
shows a zoom into the Fourier transform. The numerical solution and the analytical 
expression (|59[) exhibit besides the main peaks at u>iq and UJ20 two additional peaks, 
corresponding to the frequencies W24 (between the two main peaks) and LO23 (on the 
left of the first main peak). The two dips come from W13 and W14. 




Figure 8. Dynamics of P(t) and its Fourier transform F(oj) for positive detuning 
(£1 > Ao), for e = and £1 = 1.5Ao. The peak at ojiq dominates. No additional 
peaks are found. A very good agreement between all approaches discussed in the text 
is found. Remaining parameters are as in figure [H 

vanished. Besides, now the frequency uio is dominating the dynamics. This behaviour, 
namely that for negative detuning uj 2 q and for positive detuning co>i is dominating, was 
already found in [27] . 

We will briefly explain how one can explain this observation physically. For this we look 
at figure [9j For a detuned system (fl ^ A b ) the qubit-HO eigenstates are not symmetric 
superpositions of the states |jg) and |je). They rather asymptotically approach the 
eigenstates of the uncoupled qubit-HO Hamiltonian. In figure [1] we see that for a 
negatively detuned system (line a) the qubit-HO eigenstate \2j + 1) approaches the 
state \(j + l)g), whereas the main contribution to the state \2j + 2) will come from the 



Dissipative dynamics of a biased qubit coupled to a harmonic oscillator 



22 






Figure 9. Schematic energy spectrum for three different situations. From left to right: 
negative detuning (fi < At,), resonant case (fi = Ab) and positive detuning (f2 > Ab). 
The dashed lines show the energy levels for the uncoupled qubit-HO system (g = 0). 
The solid lines depict the eigenstates obtained by Van-Vleck perturbation theory. 



state \je). From the left diagram in figure [9] we see that the state |2) is energetically 
higher than the state |1). However, due to the Boltzmann distributed occupation of the 
oscillator, the state |0e) will be more populated than the state |lg) and consequently also 
|2) will exhibit a larger population than |1), as the latter only feels a small contribution 
from the state |0e). Thus, transitions from |2) to the groundstate are more likely to 
occur than those from |1) to the groundstate. This explains the dominance of c^o i n 
figure [7] and figure [TTJ In this case the frequency u^o ~ ^ b and lj 10 ~ fi- As far as not 
excluded by selection rules, minor peaks from transitions to the levels lying in between 
can be also seen. 

For positive detuning (line c in figured]) \2j + 1) approaches |je), while \2j + 2) is close 
to |(j + l)g). From the right graph in figure [9] we see that the state |1), being lowest 
in energy apart from the groundstate, is now also more probable to be occupied than 
1 2). Therefore, as confirmed by figure El and figure fT2| the frequency c^io is dominating 
whereas u 2 o is represented only by a small peak in the Fourier spectrum. Furthermore, 
as there are no additional energy levels between the state |1), which is most probably to 
be populated, and the ground level, other transitions than those corresponding to Uiq 
or U20 are very unlikely to occur. In figure [T2l the dip corresponding to u 2 % appears only 
very faintly. 

5.2. The biased qubit 

We will now examine a qubit being operated at finite bias. We consider the case e > 0. 
For negative bias-offset the behaviour is analogous. Again three different situations are 
taken into account: the qubit being in resonance with the oscillator (Ab = fl), negative 
(Q < A b ) and positive (f2 > A&) detuning. 

For the resonant case (Q = Ab) depicted in figure [TOl we see a similiar behaviour as for the 
unbiased qubit. Again two frequencies, uiq and u 2 o, are dominating the dynamics. Left 
to the peak at ujiq a small dip can be found in the Fourier spectrum. This corresponds to 
the transition u; 21 . For infinite time the dynamics relaxes to an equilibrium value which 



Dissipative dynamics of a biased qubit coupled to a harmonic oscillator 23 




Figure 10. Dynamics of Pit) and its Fourier transform F{uS) for the biased qubit 
being in resonance with the oscillator. Here, e — 0.5Ao and Q — At,. The remaining 
parameters are the same as for the unbiased qubit. In the Fourier spectrum the 
frequenices u>xo and u>2o dominate. At frequency u>2i a small dip can be seen. At 
lu = the spectrum exhibits a relaxation peak. 



is nonzero in contrast to the unbiased case. This can be seen in the Fourier spectrum 
through a relaxation peak at uj = 0. The peak arises because of the term 

2(Po - Poo)^^ + 2n Poo 5(Lu) (70) 

in (1591) . The first part of this sum gives rise to the negative shift of this peak. The 
reason that for the analytical solution the peak is not as strongly shifted as for the 
numerical one is technical: in order to plot the delta function in flTOj) we gave it a finite 
width, which surpresses the negative contribution of the first term in ( !70l) . Like for the 
unbiased qubit the highest energy level playing a role for the dynamics is E%\ i.e., only 
the ground and first excited level of the oscillator are of importance. In figure [11] the 
dynamics and its Fourier transform for a negatively detuned qubit-oscillator system with 
e 7^ are shown. Like for the unbiased case detuning gives raise to longer relaxation 
times for the qubit. Also in agreement with the unbiased case is the dominance of the 
frequency co> 2 o- We see that for small t the long-time solution (1571) slightly overestimates 
the maxima of the oscillations and underestimats its minima. Furthermore, we get here 
the unphysical situation that the maximum of the third oscillation in P(t) exceeds the 
value of one. The reason for that behaviour is that, by construction, we underestimate 
with ( 1571) the relaxation at short times. As for certain paramteres the term (po — poo) in 
(!57j) can become negative, it increases too fast towards the equilibrium and gives thus 
raise to the observed deviations in the short time behaviour. On a longer timescale both 
graphs agree quite well. For the case of positive detuning (f2 > At,), which is presented 
in figure [121 the upward shift of the dynamics obtained from f[57j) and ([59]) compared 
to the numerical graph of P(t) at small times is even stronger. To visualize that it is 
not a failure of the FSA approach we show in figure [11] and figure [12] additionally the 
analytical FSA solution of (142]) and (H3l) calculated in section 14.3.11 The latter agrees 
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Figure 11. Dynamics of P(t) and its Fourier transform F(oj) for negative detuning 
(il < Ab) with f2 = 0.9 and e = 0.5. Next to the numerical solution (black solid curve) 
of the full master equation and the FSA solution (|57p and (|S"9")) (red dashed curve), also 
the improved FSA solution of and (green dotted-dashed curve) are shown. 
The dynamics is dominated by L02o- The peak at lo\q is much weaker. Like for the 
resonant, biased qubit a dip is found at W21 and a relaxation peak at uj = 0. 
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Figure 12. Dynamics and Fourier transform for positve detuning (Q > Ab) for 
il = 1.5Ao and e = 0.5Ao. Like in figure QT] three different approaches are compared. 
In all three cases the frequency loiq dominates. 



very well with the numerical solution. At long time-scales and for the Fourier spectrum 
all three approaches match with each other very well. 

To conclude this paragraph we want to mention that all the results found both for the 
unbiased and the biased qubit confirm the numerical QUAPI results in [27J. 

5. 3. Symmetrized correlation function 

So far we have always considered the qubit for certain values of e and finite or zero 
detuning. In this section, we fix the oscillator frequency at Q = A . That means that 
an unbiased qubit will be at resonance with the oscillator. Changing the bias to positive 
or negative values will always lead to negative detuning, as A b > A . Figure [131 shows 
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a density plot of the Fourier transform of the symmetrized correlation function against 
the bias of the qubit and the Fourier frequency u. We consider this correlation function 
rather than P(t), as it is symmetric in the bias e. The symmetrized correlation function 





Figure 13. Fourier transform of the symmetrized correlation function plotted versus 
Fourier frequency w and qubit bias e. In the left-hand graph S(u>) is plotted in linear 
scale, in the right-hand graph 15(^)1 in logarithmic scale. The parameters are: f2 = Ao, 
g = 0.18A , K = 0.0154 and (3 = 10(^A ) -1 < Tnc wh ite dashed, horizontal line 
indicates the oscillator frequency u = CI. The other two dashed white lines correspond 
to lo = Ab and u> = Ab — O. 



is defined as follows [2]: 



1 



S(t) = §<a z (t)a z (0) + a z (0)a z (t)> - (71) 
QiHt/hp^-iHt/K Expressed in terms of the population difference P(t) this 



where a z (t) 
becomes, 

S(t) = P s (t)+ Poo (P a (t)- Poo ), 

with P s (t) and Pa.{t) being symmetric and antisymmetric in e and P(t) = P s (t) 
The Fourier transform of S(t) is defined as 



S(cj) 



dt cos(u>t) S (t) . 



(72) 

Pa(t). 

(73) 



Considering now figure [13] we see that for any bias the spectrum is dominated by 
two frequencies, namely u w and u>20- Detuning the system co> 2 o gets more and more 
important, as we could already observe in the two previous sections for the positively 
detuned systems. Furthermore, the peaks are shifted to higher frequency values and 
at uj = the relaxation peak occurs. We want to compare these results to a circuit 
QED experiment performed by Wallraff et al [13]. There the qubit is realized by a 
Cooper pair box, which is coupled to a superconducting transmission line resonator. The 
properties of the system are determined by probing the resonator spectroscopically. The 
amplitude of a microwave probe beam transmitted through the resonator is measured 
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versus the probe frequency and the gate charge of the Cooper pair box (see figure 4 in 
[13J). Via the gate charge the qubit can be detuned in situ from the degeneracy point. 
The frequency of the resonator is chosen in such a way that it is in resonance with a 
qubit being operated at the degeneracy point. For the resonant case two dominating 
frequencies, being almost equally weighted and symmetrically positioned around the 
cavity frequency, are observed. Going away from the degeneracy point the system 
becomes detuned and the frequency of the cavity dominates. The behaviour we observe 
in figure [13] is similar. However, as we are looking at the dynamics of the qubit, it 
corresponds to a spectroscopic measurement on the TLS rather than on the oscillator. 
As explained above the two lowest excited states of the coupled TLS-HO system, namely 
|1) and |2), evolve from an almost symmetric superposition of basis states {|jg), lie)} 
at resonance (e = 0) to the states |lg) and |0e) (cf the left graph in figure [9]). For 
e = the two peaks of the Rabi splitting are observed. For e ^ 0, which means negative 
detuning in this case, the peak with the lower frequency corresponding to u w approaches 
more and more the frequency Q of the oscillator, as the state |1) becomes |lg) for large 
detuning and then fhdio ~ -Ejig) — -^|0g> — Tifl- Furthermore, the transition peak at Uiq 
gets weaker as also the occupation probability of |1) decreases. At e ~ ±0.8A the 
symmetrized correlation function vanishes at c^io and increases again for higher values 
of \e\. Here, the amplitude pxo m P(f) changes its sign. In contrast the peak at uj 2 q 
becomes stronger with the detuning and approaches more and more the qubit splitting 
energy hA h , as |2) approaches |0e) and then ftw 20 ~ -Ejoe) — -Ejog) — ^A b . Additionally, 
looking at the logarithmic plot one sees around to = 0.4Ao a peak appearing, which 
corresponds to the frequency u> 2 i and is forbidden at e = 0. For large detuning it arises 
from transitions from |0e) to |lg) and therefore has the value u 2 i ~ A b — fl. The 
amplitude of this peak is very small compared to the peaks at u w and u 20 and is not 
resolved in the experiment of Wallraff et al . 

6. Comparison with the Jaynes-Cummings model 

Van-Vleck perturbation theory enabled us to find approximately the eigenstates and 
eigenenergies of the full Hamiltonian of the TLS-HO system without performing a 
rotating-wave approximation. Using those eigenstates and eigenenergies in a Born- 
Markov master equation we could calculate the dynamics of such a system under the 
influence of an environmental bath. In the following we will show how the results change 
if we neglect counter-rotating terms in the TLS-HO Hamiltonian (fl6l) for e = 0. For 
this we rewrite the interaction part in fflBl as 

Wmt = nf Qt + = -hg(d? + B + fr-flt) _ hg{a-B + a+B% (74) 

where we identified with 7^ t and Ti.™ a rotating and counter- rotating part of 7ii n t, 
respectively, and introduced the two-level transition operators d ± = \{& x ± l &y)- 
Neglecting the counter-rotating part in Hqho leads to the Jaynes-Cummings 
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Hamiltonian 



H 



JC 



h A 

-& z + mB^B - hg(a + B + a~ B ] ) 



(75) 



This Hamiltonian can be diagonalized exactly and its eigenstate and eigenvalues can 
for example be found in [46]. In order to see the effect of not taking into account 
the counter-rotating terms we diagonalize 7ijc using Van-Vleck perturbation theory. 



Looking at the formula fo the effective Hamiltonian (12.61) in Appendix B and keeping 
in mind that we set e to zero, we see that the second order contributions in g vanish 
neglecting TC^', i.e., Wi, Wo are zero in (I2T1) . Considering further the transformation 
matrix S, equations i \2A\\ and (12.51) in Appendix B show that 5* = for Ti.^ = 0. Thus, 



with f|T7|) we find that for the Jaynes-Cummings model Ti. c s is identical to Tijc and 
therefore the eigenstates of H c s are simultaneously eigenstates of Hjc- Consequently, 
one can determine from ( 122]) - (1241) the eigenstates and eigenenergies of 7ijc- The energy 
of the groundstate |0) JC = |0g) JC is Eq C = —hA h /2. For the higher states we get 



\2j + 1) JC = cos 



a 



JC 



10' + 1 )g) + sin 



a 



JC 



|je), 



(76a) 



|2j + 2) 



JC 



a 



JC 



— sin 



10' + 1 )g) + cos 



a 



JC 



lie), 



(766) 



corresponding to the eigenenergies 



E. 



2j+l/2j+2 



h 



/. , £jc 

+ o ) Q =F « jc 

2 2 cos a 



3 -I 



h 



1 



1 



(77) 



with <5jc = Ab — ft and tana:- = 2y/j + 1 1 A| / <5jc- Comparing these eigenstates and 
eigenenergies to the ones found for 7Yqho, (123 ap . (123 b\ and (|26|) . we see that the counter- 
rotating terms yield second order corrections in g not present in the Jaynes-Cummings 
Hamiltonian. These corrections give rise to a very prominent effect concerning the 
resonance condition between TLS and HO. From 5jc we find the TLS being in resonance 
with the oscillator for Q = Ab- Considering Sj in (125]) this resonance condition is shifted 
to 



n 




1 + 20 + 1 



A 2 
, o 



A 



4 r 



A, 



l + ^(2j + l)^ 2 + 0(, s 



(78) 



This second order correction to the resonance frequency due to counter-rotating terms 
is known as Bloch-Siegert shift [45]. The eigenstates (176 al) and (17661) are always a 
superposition of two basis states of the unperturbed system. This is like for the 
eigenstates (123 al) and (12361) of the effective Hamiltonian (|T7]) . However, in order to 
find the eigenstates |n) of 'Hqho we had to apply the transformation exp(— iS) on the 
effective eigenstates so that \n) is in the end a superpostion of several states of the basis 

{b'g);b'e}}. 

To calculate the reduced density matrix of the qubit-HO system described by 7Yjc and 
taking into account the influence of an environmental oscillator bath we can again use 
the Born-Markov master equation (136]) . We just have to use (176 al) and ( 176 6j) as basis 
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states and the corresponding eigenenergies. For the population difference P(t) we apply 
ffT3l). which becomes for e = 



P(t) = ^Pnnif) + ^Pnmit) 

n>m 

= ^^2(jg|n)(je|n)p nn (t) + 2 ^ ^[(je|n)(m|jg) + (je|m)(n|jg)] 9£{p„ m (t)}. (79) 



n 7 n . m 1 

n>m 



5.1. Selection rules 

Before doing a qualitative comparison between the results obtained from the original 
Hamiltonian ?^qho and the simplified form Hjc, we want to analyse which transitions 
between the different eigenstates of Hjc yield contributions to P(t). For this it is helpful 
to rewrite (176 al) and (17661) . In the following we neglect the upper index JC denoting an 
eigenstate of Hjc- For the state \n) we have three possibilities: the first one corresponds 
to n — 0. In this case the only non- vanishing component of |0) is (Og |0). Second, n 
can be an even number, which means, expressed in terms of oscillator quanta, that 
n = 2j + 2. Then 

/ n\ In — 2 \ 

|n ev ) = -sin(a^/2)| (-J g) + cos( a2if a/2)| J e>. (80) 

And third for an odd state n = 2j ' + 1 we find 

Kd> = cos(a v /2)| (—j-J g) + sin(a v /2)| J e). (81) 

It is quite easy to see that the part p n n{t) in f!79|) vanishes for any n. That means that 
the diagonal elements of the reduced density matrix yield no contributions to P(t) for 
e = and the equilibrium value of the dynamics will be also zero. Not so easy to see 
are the combinations of n and m yielding contributions to the off-diagonal part p nm (t). 
Considering transitions to the groundstate (m = 0), we find using that (0g|0) = 1 that 

Pn cv o = 2(Oe|n e v)3£{Pn cv oW} and Pn od o = 2(Oe|n od )9ft{p nod0 (t)}. (82) 

With (|80|) and (IHT!) these elements are non-zero only if n = 1 or n = 2. For the more 
general case with n ^ m and both being different from zero we have to investigate 
products of components like (je\n) (m\jg). For p ncv m cv (t) we find that 

(je\n ev )(m ev \jg) ^0 if j = Hcv and j = ^ (83) 



and that 



so that 



(je\m ev )(n ev \jg) ^0 if j = and j = (84) 



Pn cv m cv (t) ^ if \n ev - m ev | = 2. (85) 
For the case of transitions between odd states, one gets that 

(je\n od )(m od \jg) ^0 if j = n ° d and j = m ° d + 1 (86) 
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and similarly under exchange of n od and m od , so that 

Pn od m od (t) ^ if |n od - m od \ = 2. (87) 

For transitions from an even to an odd state and vice versa we have to pay attention to 
the fact that 



n ev - 2 . m od + 1 
— - — and j = . 



(88) 



(je\n ev )(m od \jg) ^ if j 
which yields the selection rule n ev — m od = 3. Further, 

O'e|m od )(n ev |jg) ^ if j = m ° d and j = ^p, (89) 
yielding m od — n ev = 1. To sum up: p n m{t) is non-zero if one of the three following 



cases is valid: \n PV — m c 



2, |n od - m 



od 



2, n ev - m od = 3 or m od - n c 



1. The 



allowed transitions are shown in figure O One sees that transitions between almost 



|o> 



Figure 14. Possible transitions between the first six eigenstates of the Jaynes- 
Cummings Hamiltonian are indicated by the red arrows. Transitions between almost 
degenerate levels are forbidden at zero bias (e = 0). 



degenerate levels are forbidden. This behaviour we have also found in section 13.31 for 
the non-dissipative dynamics resulting from 7^qho 



6.2. Comparison of the two models 

In the following we compare the numerical solution of the Born-Markov master 
equation originating from the Jaynes-Cummings Hamiltonian with the solution using 
the eigenstates of the full Hamiltonian 7Yqho; which we found by applying Van-Vleck 
perturbation theory. The fixed parameters we use are in units of Ao: e = 0, g = 0.18, 
k = 0.0154 and (3 = 10. The oscillator frequency Q is varied. For all the three possible 
cases (positive, negative and zero detuning) one notices from figures [15] - [17] that the 
Jaynes-Cummings approach underestimates the dephasing rate Tio (means a larger peak 
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at frequency uio) and overestimates the rate r 2 o (smaller peak at u 2 o) compared to the 
approach with the full Hamiltonian. 

The case of positive detuning (Q > Aq) is shown in figure [151 As here the dynamics 
is dominated by the frequency u w , the equilibrium value is reached on a too long time 
scale using the Jaynes-Cummings Hamiltonian. On the contrary, for negative detuning 
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Figure 15. Dynamics P(t) of the population difference and its Fourier transform 
F(u>) for positive detuning (f2 = 1.5Aq). The parameters are in units of Ao: e = 0, 
g = 0.18, k = 0.0154 and (3 — 10. The red dashed line shows the solution obtained 
numerically from the master equation using the eigenstates and eigenenergies from the 
full Hamiltonian Hquo- The solid line shows the results obtained from the Jaynes- 
Cummings Hamiltonian. 



(Q < A ), which is shown in figure [16] and where u 2 o dominates, the equilibrium value 
is reached too fast within the Jaynes-Cummings approach. Furthermore, considering 
the graph of the Fourier transform we find that small contributions which come from 
higher level transitions, and which have already been discussed in section 15.11 are not 
caught be the Jaynes-Cummings approach. For the resonant case (Q = A ) we find 
that the Jaynes-Cummings method predicts Uiq to be slightly dominating whereas the 
approach starting from Hqko results in co> 2 o being dominating. The reason for this 
discrepancy is that due to the counter-rotating terms we have for 7Yqho no symmetric 
or antisymmetric superposition of the unperturbed eigenstates at Q = Aq in contrast 
to the Jaynes-Cummings model. 

To conclude this section we can say that for an unbiased TLS-HO system the Jaynes- 
Cummings model gives a good insight in the qualitative behaviour of P(t) both for 
a slightly detuned and a non-detuned system. However, it under- or overestimates 
dephasing times for the system. Furthermore, we find taking into account counter- 
rotating terms in H.quo that at Q = A^ the dressed eigenstates are not a symmetric or 
antisymmetric superposition of the uncoupled states. Moreover, the effects of transitions 
between states of different manifolds are neglected. 
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Figure 16. Dynamics P(t) of the population difference and its Fourier transform F(ui) 
for negative detuning (f2 = 0.75Ao). The remaining parameters are the same as in 
figure[T5] The red dashed line shows the solution obtained numerically from the master 
equation using the eigenstates and eigenenergies from the full Hamiltonian 7Yqho • The 
solid line shows the results obtained from the Jaynes-Cummings Hamiltonian. 




Figure 17. Dynamics P(t) of the population difference and its Fourier transform F(u) 
for the resonant case (fi — Ao). The remaining parameters are the same as in figure [151 
The red dashed line shows the solution obtained numerically from the master equation 
using the eigenstates and eigenenergies from the full Hamiltonian 7Yqho- The solid 
line shows the results obtained from the Jaynes-Cummings Hamiltonian. 



7. Conclusions 



In conclusion, we discussed the dynamics of a biased and unbiased TLS coupled through 
a harmonic oscillator to an environmental bath described by an Ohmic spectral density. 
In particular, we examined the regime of weak damping and moderate coupling between 
oscillator and TLS. An equivalent description of our system is provided by the spin- 
boson model with a structured spectral density. In contrast to many other works in 
this field, our starting point was not the Jaynes-Cummings Hamiltonian, but a more 
general one given in (JT]), where no initial rotating wave approximation has been applied. 
In section 12.21 we provided with ( Tl3|) a formal expression for the population difference 
and showed in section I3~T1 how the Hamiltonian of the coupled qubit-HO system can be 
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diagonalized approximately using Van-Vleck perturbation theory. This approach is valid 
both for the oscillator being in resonance with the qubit (Q = At,) and for finite detuning 
(\S\ 7^ 0). In section [32] an analytical expression for the non-dissipative dynamics was 
provided up to second order in the qubit-HO coupling g taking into account the infinite 
Hilbert space of the system. For low temperatures (k&T < hQ, hA\,) we truncated the 
Hilbert space and found that transition processes between the groundstate and the two 
first excited energy levels of the qubit-HO system dominate the dynamics (section 13. 3j) . 
In section HJthe influence of the bath was taken into account by solving the Born-Markov 
master equation for the density matrix of the qubit-HO system. To do this analytically 
we considered two variants of the secular approximation: first, in section 14.3.11 the full 
secular approximation (FSA), where all fast oscillating terms are neglected, and second, 
in section 14.3.31 the partial secular approximation, where attention was paid to the 
fact that the first two excited energy levels are almost degenerate. Using an ansatz for 
the long time dynamics in section I4.3.2[ we could provide a general expression for the 
relaxation and dephasing rates of the qubit, showing that the relaxation time can be 
enhanced by detuning the oscillator into the off-resonant regime. It was found that all 
three approaches agree quite well with the numerical solution of the master equation. 
The dynamics of both a biased and an unbiased qubit were intesively studied for zero and 
finite detuning in section [51 The results agree qualitatively with the numerical findings 
within the ab-initio QUAPI approach [27]. Furthermore, in section I57T1 a good agreement 
with the results of the weak damping approximation performed in [31] for a symmetric 
spin-boson model was found. Besides, we saw that at resonance (ft = At,) the first two 
excited qubit oscillator states are not a symmetric or antisymmetric superposition of 
the states |0e) and |lg), as predicted by the Jaynes-Cummings model, and thus could 
give an explanation for the differently weighted peaks of uj 10 and cg> 20 i n the Fourier 
spectrum of the dynamics. We further could explain the dominance of frequency u^o 
in the case of negative detuning (Q < At,) and of frequency uio for positive detuning, 
respectively. Moreover, we showed that for large negative detuning co> 2 o approaches the 
energy splitting At, of the qubit, whereas uj\o approximates the oscillator frequency 
Q. This behaviour agrees nicely with spectroscopic experiments performed on a circuit 
QED architecture [13] . In section [6] we compared our results for an unbiased system to 
the ones obtained starting with the Jaynes-Cummings model. We visualize the effects 
of the counter-rotating terms in the Bloch-Siegert shift of the resonance frequency and 
in contributions of states with larger oscillator number to the TLS dynamics. Apart 
from this the Jaynes-Cummings model and our approach agree quite well for the non- 
dissipative case. Also for the dissipative case at resonance an initial RWA represents 
a good approximation to our starting Hamiltonian and seems to be favourable as it is 
analytically exactly diagonalizable. For detuned systems, however, we find discrepancies 
concerning relaxation and dephasing times and a RWA becomes less appropriate to give 
precise results. Thus, we think that our approach represents an improvement as it is 
valid in a wider parameter range avoiding an initial rotating wave approximation. To 
our knowledge it provides for the first time analytical results for the dynamics of an 
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unbiased and biased qubit coupled to a structured environment being valid both in 
the resonant and off-resonant regime. Furthermore, due to the generality of the qubit- 
oscillator model, we expect our results to be of interest for a wide range of experimental 
applications. 
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Appendix A. Formula for the dynamics 

Here, we show how to derive ([TBI given in section 12.21 In order to trace out the 
oscillator degrees of freedom we transform p nm (t) into the basis {|jg); |je)} and by 
using p nm = p* mn we find 

PjgjM = (jg|p(*)Ug) = ^2(jg\n) 2 pnn{t) + ^ ^{Pnrn{t)}(]g\n)(m\jg) (Lid) 

n ™. m 

Pje,je(t) = (je\p(t)\je) = ^2(je\n) 2 p nn (t) + ^ ^{pnm(t)}(je\n}(m\je}. (1.1b) 

n n,m 

Performing the trace over the oscillator 

oo oo 

Pred;gg(*) = (g|Pred(*)|g) = ^^ij g\ n } 2 Pnn(t) + ^ ^{Pnm(t)} (j g\n) (m\j g) , (1.2a) 

j=0 n j=0 n - m 

oo oo 

Prcd (t) = (e\ Pied (t)\e) = ^5Z(je|n) 2 P nn(t) + ^{Pn m (t)}(je\n)(m\je). (1.2b) 

j=0 n j=0 ». m 

Similarily, we find for the off-diagonal elements of the reduced density matrix 

oo 

Pred;eg(t) = ^^{j e\n) p nn (t) (n\j g) 
j=0 n 

^ oo 

+ 2 S Yl [(j e \ n )Pnm(t)( m \j&) + (j e \ m )P*nm(t)(jg\n)} , (l-3o) 

j — Q n,m 

oo 

Pred;ge 

(*) = J2J2(jg\ n )Pnn(t)(n\je) 

j=0 n 

+ 2EE [(jg\ n )Pnm(t)(m\je) + (jg\m)p* nm (t)(je\n)] . (1.3b) 

j — Q n,m 
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Using (jSJ) and flSD in (TTT|) we can express P(£) in the energy basis, yielding 

P(t) = COS 6 [p r ed;gg(*) ~ Pred;cc(£)] + Sm 6 [p re d;ge(^) + Pred;eg(£)] 

oo 

= cos(6) ^ [(jg\n) 2 - (je\n) 2 } p nn (t) 



j=0 n 

oo 



+ [(js\n)(m\jg) - (jg\n}(m\jg)} 5£{p nm (£)} 

j = Q rt,m 

oo 

+ sin(0) (2 ^2(Mn)Pnn{t)(n\jg) 



j=Q n 



+ Y1 2 [0' e l n )( m b'g) + (j e |w)(jg|n)]3£{p nm (t)} 



j=0 n,m 



n . in 

n>m 



+ sin 9 



:i-4) 



With (J28|) we can write 

-P(^) = ^Pnn(^) + ^ Pnm(*) COSU; nm t, (1.5) 



where 

Pnn(^) = |cos9 (jg\n) 2 - (je\n) 2 + 2sme(jg\n) (je\n)jp nn (t), (1.6a) 

3 

p nm (t) = 2^|cos6 (jg\n)(m\jg) - (je\n)(m\je) 



(je\n)(m\jg) + (je\m)(n\jg)^ }5?{p nm (£)}. (1.6b) 



Appendix B. Van-Vleck perturbation theory 

Let us consider the Hamiltonian 

H = H + V (2.1) 

consisting of the free Hamiltonian 7io and a small perturbation V ~ g, which is 
proportional to the coupling constant g. Additionally we assume, that the energy levels 
Ej t0l of Tio are grouped into manifolds, with a being the index of the manifold and i 
is used to distinguish between different energy levels within the same manifold. The 
energy levels E^ a are eigenenergies of TCo: 

H \j,a) = E Jta \j,a). (2.2) 

Through the perturbation V different manifolds are coupled together. As long as the 
coupling g is small, namely that | (j, a\V\j, j3) | <C \Ej >a — Ej t p\ for a ^ j3, also the energy 
levels of the total Hamiltonian 7i are clustered into manifolds. Using the transformation 
7~(-cE — e l5 7^e~ 1,s , we construct an effective Hamiltonian 7i c s, which acts only within the 
individual manifolds; i.e., (j, a\TC e s\j, ft) = for a ^ (3, and has the same eigenvalues 
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as 7i within the manifolds. We expand S and 7i e ff in terms of the small parameter g up 
to second order: 

S = s w + s m + 0{g 3 ) and HcS = n w + n w + n m + 0{g 3 ) 



(2.3) 



For calculating S 1 * 1 / 2 ) and H.*J we use that (J, a\H e Jf \j, (3) = for a ^ (3 
and furthermore choose that S has no matrix elements within a manifold, namely 
(j,a\iSW\j,a) = 0. Now, one can iteratively calculate S and 7i e fr order by order. 
For details see [JT]. Here, we give only the results. For the transformation one has 

(j,a\V\j,i3) 



{j,a\iSM\jJ) 



for a 7^ P, 



and 



aai^ij^> = - Yl 

1 



(j,a|\/|fc,7)(fc,7|^|j,/3) 



+ 



(j,a\V\k,P)(k,P\V\j,P) 



Ek,f3 — Ej^ a 



{j,a\V\k,a)(k,a\V\j,fl 



+£ 

k Ej,P ~ Ej, a E kyCt - Ej£ 

The effective Hamiltonian is up to second order 
(i, a\H e ff\j, a) = E jta 5ij + (i, a\V\j, a) 
+]- (i,a\V\k,j)(k,j\V\j,a) 



for a 7^ j3. 



k,~f=£a 



+ 



1 



Ei y a E k ^i Ej^a Efory 



+ 0(g 3 



In the case of the Hamiltonian 7"^qho the first order matrix elements are 



59(1) 

w gjgj+l 



6J e 3 + l V J 



g, 



A b tt 
A 



and for the second order contributions 



;9(2) 

e i e i+2 



i»5p .„ . 



]C(2) 

gjgj+2 



]C(2) 

gj e J+2 



2v/(j + 1)(j + 2) 



eA r 



A 2 fi(2ft - A 6 



V(i + l)(j + 2) 



An 



n gA 
2 l J j A 2 fi(A b + 0) 



2A 2 fi(A 6 + fi) 

2 



9 



9 



A 2 



iva+i)a+2) AW+A6) 



A 2 b Q{A b + Q){A b + 2Q)- 



(2-4) 



;2.s) 



(2.6) 



(2.7a) 
(2.76) 

(2.7c) 

(2.8a) 
(2.86) 
(2.8c) 
(2.8d) 
(2.8e) 
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where e. g. 

Sif_ iej = {(j-l)e\S\je). (2.9) 

These matrix elements change sign under index transposition and all other matrix 
elements vanish. Finally, we get the transformation up to second order in g: 

e ±iS = 1 ± iS^ ± iS® + Us^iS^ + 0(g 3 ). (2.10) 

Appendix C. Oscillator matrix elements 

Here, we give the matrix elements X nm specified in section 14.21 
X 2 j+i,2j+i = -2L Q cosa>j + ^/j~+TL~ osc sma j , 
X 2 j+i,2j+2 = 2L sinaj + \J j + lL~ osc cosa,-, 

X 2 j+i,2j+3 = L q cos(oj/2) sin(a i+ i/2) + y/j + 2(1 + L osc ) cos(a.,/2) cos(a j+ i/2) 

+Vj + 1(1 - L osc ) sin(o;j/2) sin(a j+ i/2) 
X 2 j+i,2j+4 = L q cos(a i /2) cos(a j+ i/2) - y/ j + 2(1 + L osc ) cos(a i /2) sin(a i+1 /2) 

+ + 1(1 - Lose) sm(aj/2) cos(a i+ i/2) 
X 2j+ i,2j+5 = Vl + 2L cto SC cos(ay/2) sin(a i+2 /2), 
X 2j+ i,2j+e = \fj + 2L+ osc cos(a j /2) cos(a i+2 /2), 
X 2 j+2,2j+2 = 2L cos aij - \/j + lLq OSC sinaj, 

X 2 j+2,2j+3 = -L q sin(aj/2) sin(a j+ i/2) - a/j + 2(1 + L osc ) sin(<x,-/2) cos(a i+ i/2) 

+ \Jj + 1(1 - L osc ) cos(a i /2) sin(a i+ i/2) 
X 2 j+2,2j+4 = -L q sm(aj/2) cos(a j+ i/2) + yj j + 2(1 + L osc ) sin(a,/2) sin(a j+ i/2) 

+ a/j + 1(1 - L osc ) cos(a,/2) cos(a j+ i/2) 
X 2 j+2,2j+5 = -y/j + 2L tosc sin(aj/2) sin(a j+2 /2), 
X 2j +2,2j+6 = - v / 7T2L+ osc sin(a j /2)cos(a i+2 /2). 

Matrix elements including the groundstate are given separately because of the special 
shape of |0): 

Xo,o = — 2L , 

X 0) i = sin(a /2)L q + cos(ot /2)(l + L osc ), 
X ,2 = cos(a /2)L q - sin(a /2)(l + L osc ), 
X ,3 = sin(a 1 /2)L+ osc , 
X 0A = cos(a!/2)L+ osc . 

All other matrix elements are zero. 
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Appendix D. Rate coefficients for the off-diagonal density matrix elements 

Here, we give the rate coefficients of the master equation for the reduced density matrix 
elements. They are: 

Ak 1 

£01,01 = ^^o( 2cosa o - cos 2 « - 1) - 2^00,11, (4.1) 

Ak 1 

£02,02 = -^jLl(-2 cos a - cos 2 a - 1) - -£00,22, (4.2) 

£03,03 = ]^j L o( 2 cos «i - cos2 - 1) - ^£11,33 - ^£22,33, (4.3) 
Ak 11 

£04,04 = ^p L o(~ 2 COS "1 ~ COs2 «1 - 1) - -^C-UM ~ T^ 22 ' 44 ' ( 4 - 4 ) 

16/t 1 1 

£12,12 = — ftjJ L o cos2 a ° ~~ 2^ 00,11 ~~ 2^ 00 ' 22 ' 

£13,13 = -^L 2 (cosa - cosoi) 2 - ^£00,11 - ^£11,33 - ^£22,33, (4.6) 

£14,14 = -j^Ll(cosa + cosoi) 2 - ^£00,11 - ^£11,44 - ^£22,44, (4.7) 

£23,23 = -^Ll(cosa + cosoi) 2 - ^£00,22 - ^£11,33 - ^£22,33, (4.8) 

£24,24 = -j^jLl(cosa - cos«i) 2 - ^£00,22 - ^£n,44 - ^£22,44, (4.9) 

£12,12 = —^t^LI cos 2 «i — -£11,33 — -£22,33 — -£11,44 — -£22,44, (4-10) 
Ak 

£oi,02 — ^(^00^12 — -^12^22) — G {u)q 2 )Nq 2 XqiXq 2 — G (uJ 12 ) Ni 2 XuXi 2 

-G{uo 32 )N 32 X 13 X 23 - G(u 42 )N 42 X 14 X 24 , (4.11) 

Ak 

£02,01 — ^(^00^12 — -^12^11) ~~ G{luqi)Nq 1 XqiXq2 — G{lo 2 i)N 2 iX 22 Xi 2 

-G(u 31 )N 31 X 13 X 23 - G(uj 41 )N 41 X 14 X 24 , (4.12) 

Ak 

£13,23 — ^(^33-^12 _ Xi 2 X 22 ) — G(loi 2 )N 12 (XuXi 2 — X 12 X 33 ) 

-6^02)^02*01*02 - G{uj 32 )N 32 X 13 X 23 - G(u 42 )N 42 X 14 X 24 , (4.13) 
Ak 

£23,13 — ^(*33*12 — *12*ll) — C(c<J2l)iV"2l(X22^12 — *12*33) 

-G(uj 01 )N 01 X 01 X 02 - G(u 31 )N 31 X 13 X 23 - G(uj 41 )N 41 X 14 X 24 , (4.14) 
Ak 

£14,24 = ^(*44*i2 — Xi 2 X 22 ) — G(u>i 2 ) N 12 (XuX 12 — X 12 X 44 ) 

—G(u>o 2 )N Q2 X 01 X 02 — G(u 32 )N 32 X 13 X 23 — G(lo 42 )N 42 X 14 X 24 , (4-15) 
Ak 

£24,14 = jrjj(X 44 Xi 2 — Xi 2 Xu) — G(u 2 i) N 2 i(X 22 Xi 2 — X12X33) 

-G(u;oi)*oi*oi*02 - G(u 31 )N 31 X 13 X 23 - G(u 41 )N 41 X 14 X 24 . (4.16) 
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Appendix E. Diagonal reduced density matrix elements 



The solutions of the FSA master equation (1511 for the diagonal elements reads: 

^00W = <%) + °11 + °22 + ^33 + °44 

_ p -7r£ 00 ,iit ( , _0 Al.33 n . £ H,44 

c i u n -r ^33 

/* i r \ r ' 44 



~£oO,ll + £ll,33 + £22,33 — £oO,ll + £11,44 + £22,44 
^C 00}22 t(n , n ^2,33 ^22,44 



\ °22 ' °33 , , n ' °44 



_|_g— 7r(£ll,33+-C22,33)t ( j0^ / 



£00,22 + £ll,33 + £22,33 — £00,22 + £ll,44 + £22,44 

£00,22 — £ll,33 £ll,33 



— £00,22 + £ll,33 + £22,33 — £00,11 + £ll,33 + £22, 



33 

I e -n(Cn M +c 2 2 M )t a a ( £00,22 - £n,44 + £n,44 1 

\— £00,22 + £11,44 + £22,44 — £00,11 + £11,44 + £22,44- 

„ (+\ — o-kCoo, lit 

_ -7r(£oO,ll+£ll,33+£22,33)t n -0 £ll,33 



a 



33 



"£00,11 + £11,33 + £22, 



33 



— 7r(£oo, 11 44+.C22, 44 )* n -0 



£ 



. r , (5-2) 



— £oO,ll + £ll,44 + £22, 



11 



22 

r 

-7r(£oO,22+-Cll, 33+^22, 33)*^-0 22,33 



33 



"£oO,22 + £ll,33 + £22, 



33 



(3 -7r(£ 0,22+£ll,44+£22,44)t_0 £22,44 /- „n 

— ^00,22 1 Ml,44 + ^22,44 

^(t) = e — ( £ ",33+£22,33)V0 3) (5.4) 

a 44 (t) = e-^^+^.^V^. (5.5) 
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